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1.0  Summary 


This  report  describes  the  March  1994  Arctic  deployment  undertaken  by  the  Acoustic 
Telemetry  Group  of  WHOL  The  deployment  was  a  part  of  the  1994  Sea  Ice  Mechanics 
Initiative  (SIMI)  project  and  was  based  at  the  west  SIMI  camp,  approximately  150  Nau¬ 
tical  miles  north-east  of  Prudhoe  Bay,  Alaska.  The  goal  of  the  deployment  was  to  install 
a  network  of  six  high-performance  acoustic  modems,  developed  at  WHOI,  and  to  obtain 
a  data  set  demonstrating  the  communications  and  acoustic  monitoring  capabilities  of  the 
network. 

The  six  modems  in  the  network  were  deployed  over  an  area  of  22  square  km  and  com¬ 
municated  via  radio  Ethernet  with  a  computer  at  the  SMI  camp.  Each  modem  had  a 
global  positioning  system,  an  acoustic  source  and  an  8  element  receiving  array.  The  net¬ 
work  was  operated  in  a  round-robin  broadcast  mode  (i.e.,  each  modem  in  turn  transmit¬ 
ted  a  packet  of  data  while  the  others  received).  The  transmissions  were  5000  bits-per- 
second  QPSK  with  a  15kHz  carrier.  An  extensive  data  set  including  raw  acoustic  data, 
source  localization  information,  and  modem  position  was  collected  during  the  deploy¬ 
ment. 

An  additional  function  of  the  acoustic  network  was  to  communicate  with,  and  track,  the 
Odyssey,  an  autonomous  imderwater  vehicle  operated  by  the  MIT  group  at  the  SMI 
camp.  To  this  end,  the  Odyssey  was  equipped  with  a  Datasonics  modem  configured  for 
periodic  QPSK  transmission  to  the  network.  A  data  set  was  obtained  from  which  both 
the  up-link  communication  and  localization  capabilities  of  the  network  can  be  deter¬ 
mined. 
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2.0  introduction 


The  Arctic  deployment  took  place  from  March  21  to  March  30,  1994,  at  the  west  Sea 
Ice  Mechanics  Initiative  (SIMI)  camp  which  was  situated  approximately  150  Nautical 
miles  north-east  of  Prudhoe  Bay,  Alaska.  The  WHOI  Advanced  Engineering  Laboratory 
(ART.)  personnel  involved  in  the  deployment  were  Josko  Catipovic,  Ed  Denton,  David 
Herold,  and  Mark  Johnson. 


The  aim  of  the  deployment  was  to  install  and  operate  a  network  of  six  high-performance 
acoustic  modems  developed  by  AEL.  As  shown  conceptually  in  Fig.  1,  each  modem 
was  equipped  with  an  acoustic  source,  an  acoustic  receiving  array,  and  a  radio  Ethernet 
link  to  a  Sun  SPARC  computer  at  the  base  camp.  The  base  camp  computer  was  used 
both  to  collect  data  from  the  modems  and  to  control  the  modem  functions. 


BASE  CAMP 


10m  Mast 


6  Modems  remote 
from  base  camp 


Fig.  1  Conceptual  diagram  of  the  Arctic  deployment. 


The  modem  hardware  shown  in  Fig.  2,  comprises  a  personal-computer  (PC)  system 
controller,  two  high-speed  digital  signal  processors  (DSPs),  analog-to-digital  converters 
(ADCs),  digital-to-analog  converters  (DACs),  and  a  global  positioning  system  (GPS) 
together  with  a  number  of  diagnostic  peripherals.  The  unit  is  battery-powered  and  incor¬ 
porates  both  power  and  temperature  sensing  and  control  circuitry.  Further  details  of  the 
modem  hardware  and  software  may  be  found  in.  From  the  viewpoint  of  the  network 
operator,  each  modem  is  associated  with  a  DOS  window  on  the  base  camp  computer, 
from  which  all  aspects  of  the  modem  can  be  directly  controlled. 

The  modem  network  deployed  in  the  Arctic  was  configured  to  operate  in  a  broadcast 
mode  in  which  one  modem  would  transmit  and  the  remaining  modems  in  the  network 
would  receive  the  transmitted  data.  The  transmissions  used  quadrature  phase  shift,  key¬ 
ing  (QPSK)  modulation  with  a  earner  frequency  of  15kHz  and  a  data  rate  of  5000bits- 
per-second  (bps).  Following  a  transmission,  each  modem  logged  the  received  signal 
(demodulated  and  decimated  to  reduce  storage  requirements)  and  its  precise  arrival  time 
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Fig.  2  System  Diagram  of  Acoustic  Modem  with  associated  peripherals 

together  with  diagnostic  data  regarding  the  power  consumption,  temperature,  and  GPS 
position  of  the  modem.  A  report  was  also  available  from  each  modem  indicating  the 
performance  of  the  receiver  algorithm  in  decoding  the  transmitted  signal.  Experimental 
results  were  collected  by  the  base  camp  computer  and  displayed  using  Matlab.  The 
parameters  used  in  each  modem  for  transmission  and  reception  (e.g.  source  level,  syn¬ 
chronization  codes,  equalizer  settings,  etc.),  could  be  readily  accessed  and  changed 
from  the  base  camp  computer. 

An  additional  function  of  the  Arctic  network  was  to  demonstrate  communication  with, 
and  trackmg  of,  the  Odyssey,  an  autonomous  underwater  vehicle  ( AUV)  operated  by  the 
MU  group  at  the  SMI  camp.  The  Datasonics  modem  used  on  the  Odyssey  was  config¬ 
ured  solely  for  periodic  transmission  preventing  two-way  communication  with  the 
AUV.  However,  data  obtained  with  the  network  acting  as  a  receiver  demonstrates  both 
the  up-link  communication  and  localization  capabilities  of  the  network. 

In  the  following  section,  some  aspects  of  the  design  of  the  communication  network  are 
discussed.  A  summary  of  the  data-logging  experiments  performed  with  the  network  is 
then  given  in  Section  4  followed  by  a  description  of  the  logged  data  formats  and  struc¬ 
tures  in  Section  5.  Complete  experimental  descriptions  are  given  in  Section  6  and  Mat- 
lab  and  C  programs  used  to  extract  information  from  die  acquired  data  base  are 
provided  in  Section  7.  The  report  concludes  with  references  and  appendices. 
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3.0  Communication  System  Description _ 

In  the  following  sub-sections  details  of  the  transmitted  signals  and  the  transmit  and 
receive  transducers  are  given. 

3.1  Acoustic  Transducers 

The  vertical  receive  array  used  with  each  modem  comprised  8  omni-directional  trans¬ 
ducers  each  separated  by  Im.  The  transducers  were  numbered  1-8  with  element  1  being 
the  lowest.  By  convention,  element  1  was  used  for  single  channel  reception  while  ele¬ 
ments  1,3,5,  and  7  were  used  for  four  channel  reception  (i.e.,  an  effective  inter-element 
spacing  of  2m). 

The  transmit  sources  were  Datasonics  Ar-18DT  15kHz  directional  transducers  with  a 
beamwidth  of  approximately  +/- 15  degrees  from  the  horizontal.  The  usable  frequency 
band  of  these  transducers  is  between  8kHz  and  22kHz. 

Modem  sites  were  chosen  with  low  ice  thickness  (i.e.,  <3m)  and  the  sources  and  hydro¬ 
phone  arrays  were  lowered  through  holes  in  the  ice  to  a  depth  of  about  lOQm. 

3.2  Transmitted  Signal  Description 

All  transmissions  were  5000bps  QPSK  centered  at  15.151kHz.  The  transmitted  packet 
consists  of  three  sections:  first,  a  13  chip  Barker  code  sequence  was  transmitted  at  a  rate 
of  1  chip  per  symbol  to  enable  synchronization  at  the  receiver.  This  was  followed  by  a 
period  of  silence  (typically  of  length  20  symbols),  required  to  enhance  receiver  Barker 
detection  in  multipath  environments.  The  third  section  of  the  transmitted  packet  con¬ 
tained  the  QPSK-modulated  data.  The  number  of  data  symbols  in  a  packet  was  typically 
1(X)0, 2000,  or  1(X)00  depending  upon  the  number  of  receiver  channels  used. 

The  settings  used  by  the  transmitter  are  contained  in  the  txm.m  file,  a  text  file  used 
directly  by  the  DSPs  in  each  modem.  An  example  txm.m  file  is  given  below: 

%  txm.m  for  10000  symbol  QPSK  packet  transmission 

GPS_EXIST=1  %  Enables  GPS-based  timing 

time(l)  %  Starts  GPS  pulse-per-second  timer 

file  =  Ttrain.bin*  %  Binary  data  file  to  be  transmitted 

map  =  [i4j;-l-j;-l+j;l-j];  %  Symbol  mapping  for  QPSK 

code=  [-1;1;-1;1;-1;-1;1;1;-1;-1;-1;-1;-1] ;%  13  chip  Barker  code 
txf=tl;l;l;l;l;l;l;l;l;l;l;l]: 
nullf  =  [0;0;0;0;0;0;0;0;0;0;0;0]; 

txf  =  [txf;txf;nullf;nullfl ;  %  Passband  transmit  filter 

rate  =  [2500;  15 151  ;60000;0.5;0;20;  10000:200;  1.5 ;0.0];%  Options  for  transmitter 

transmit(file,map,code,txf.rate)  %  Transmit  a  padcet 

quit  %  Return  to  DOS  shell 


4of50 


Report  on  the  Acoustic  Network  Arctic  Deployment,  March  1994 


Communication  System  Description 


The  numbers  in  the  rate  vector  defined  above  control  a  number  of  aspects  of  transmitter 
performance  as  follows:. 

Option:  Units  Place  in  rate  vector 


data-rate 

symbols/sec 

1 

carrier  frequency 

Hz 

2 

samphng-rate 

Hz 

3 

amplitude 

- 

4 

output  channel 

- 

5 

silence  period  length 

symbols 

6 

packet  length 

symbols 

7 

block  size  (internal) 

symbols 

8 

Barker  code  magni¬ 
tude  (real  and  imag.) 

- 

9.10 

NOTE:  No  pulse  shaping  or  transmit  filtering  was  applied. 

3.3  Receiver  Description 

The  settings  used  by  the  receiver  are  contained  in  the  rxdLm  file,  a  text  file  used  directly 
by  the  DSPs  in  each  modem.  An  example  rxdLm  file  is  given  below: 

%  rxdl.m  for  1000  symbol  QPSK  packet  recq)tion 

file  =  ’out.da’  %  data  file 

map  =  [l+j;-l-j;-l+j;l*j];  %  qpsk  constellation 

code 

;  %  barker  sequence 

rxf  =  [-0.0013;-0.0011;-0.0008;0,0;0.0019;0.0053;0.0107] ; 
rxf  =  [rxf;0.0182;0.0277;0.0388;0.0508;0.0626;0.0734;0.0819]; 
rxf  =  [rxf;0.0874;0.0893;0.0874;0.0819;0.0734;0.0626;0.0508] ; 
rxf=[rxf;0.0388;0.0277;0.0182;0.0107;0.0053;0.0019;0.0]; 
rxf  =  [rxf;-0.0008;-0.0011;-0.0013] ;  %  receive  filter  settings 

% 

N  =  [30,1:30.0] ;  %  equalizer  settings 

pli  =  [0,0011.-0.00L0;L-2.1] ;  %  pll  settings 

rate  =  [2500;15151;60000;35;1000;0.99;LO;2;0;0.001;0.05;100;10;20;350;50.0]  ; 

‘  receive(file,map,code,rxf,'qtrain.dat',N,pll,rate);  %  receive  a  packet 
quit  %  exit 

The  key  features  of  this  file  are:  1)  the  receive  filter  (rxf)  which  is  a  3 1  tap  linear-phase 
windowed  filter,  2)  N  which  tells  the  number  of  feedforward  and  feedback  parameters  in 
the  equalizer  as  well  as  the  channel  number  the  data  was  received  on  (N(ff  parameters, 
channel;  ff  parameters  .channel; ...  ;fb  parameters,  fb  channel)).  The  critical  components 
in  the  rate  vector  are:  2500  symbol  rate,  15151  carrier  frequency,  60000  sample  fre¬ 
quency,  1000  packet  length  in  symbols,  and  2  which  is  the  inverse  of  the  fractional  spac¬ 
ing  and  determines  the  data  rate  of  the  logged  sampled  data  (i.e.  2  x  2500  =  5(X)0 
samples  per  second).  The  xtual  symboling  rate  is  1/24  of  the  sample  rate  which  is 
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determined  by  2e6Pi,  where  X  is  an  integer.  In  this  case  X  =  33  and  the  sampling  rate  is 
60606  Hz.  thus  the  actual  symbol  rate  is  60606/24  =  2525  symbols  per  second. 


4.0  Experiment  Summary 


4.1  Summary  of  Data  Capture  Experiments 

A  total  of  9  data  capture  experiruents  were  performed  between  March  26  and  March  30 
yielding  about  SOOMbytes  of  data.  The  experiments  are  summarized  below  and 
described  in  more  detail  in  Section  6.  The  experiment  code  given  for  each  experiment 
provides  a  cross-reference  with  the  log  book  of  the  Arctic  deployment.  All  times  given 
in  this  report  are  local  (central  Alaskan)  time  which  is  9  hours  behind  GMT. 

Expt.  1:  March  26, 8:10  Round-robin  transmissions  in  a  four  node  network;  lOK 

symbol  packets  (also  some  decoded  IK  symbol  packets);  single  channel  recep¬ 
tion.  This  data  set  was  taken  the  morning  after  a  substantial  lead  opened 
between  Sites  3.2  and  Sites  1.6.  Experiment  code:  rr, 

Expt.  2:  March  21, 14:03  Transmissions  from  a  Datasonics  modem  on  the  MTT 

Odyssey  received  by  a  five  node  network;  IK  symbol  packets;  single  channel 
reception.  Experiment  code:  mitL 

Expt.  3:  March  27, 22:38  Roimd-robin  transmissions  in  a  four  node  network  with  a 
change  of  transmitter  every  hour;  5K  symbol  packets;  single  channel  recep¬ 
tion.  Experiment  code:  on27_3, 

Expt.  4:  March  28, 10:12  Roimd-robin  transmissions  in  a  three  node  network  with  a 
change  of  transmitter  every  hour;  lOK  symbol  packets;  single  channel  recep¬ 
tion.  The  Datasonics  modem  on  the  MTT  Odyssey  was  transmitting  during  the 
final  hour  of  this  experiment.  Experiment  code:  on28_3. 

Expt.  5:  March  28, 22:00  Transmissions  from  a  Datasonics  modem  on  the  MIT 
Odyssey  received  by  a  six  node  network;  IK  symbol  packets;  four  channel 
reception.  Experiment  code:  mit2, 

Expt.  6:  March  28, 22:22  Roimd-robin  transmissions  in  a  six  node  network  with  a 
change  of  transmitter  every  hour;  2K  symbol  packets;  four  channel  reception. 
Experiment  code:  on28_3a, 

Expt.  7:  March  29,  9:57  Round-robin  transmissions  in  a  six  node  network  with  a 
change  of  transmitter  every  hour;  2K  symbol  packets;  four  channel  reception. 
Experiment  code:  on29_3, 

Expt.  8:  March  29, 22:40  Round-robin  transmissions  in  a  six  node  network  with  a 
change  of  transmitter  every  hour;  lOK  symbol  packets;  single  channel  recep¬ 
tion.  Experiment  code:  on29_3a. 
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Expt.  9:  March  30, 10:05  Simultaneous  transmissions  from  two  modems  to  remain¬ 
ing  modems  in  a  six  node  network;  25K  symbol  packets;  four  channel  recep¬ 
tion.  Experiment  code:  collide. 

4.2  Modem  Site  and  Package  Naming  Convention 

Six  modem  sites  were  selected  for  the  network  nodes  covering  an  area  of  approximately 
22  square  km  roughly  centered  on  the  base  camp.  Acoustic  transducers  were  perma¬ 
nently  installed  at  each  site.  However,  the  modem  electronics  packages  (contained  in 
insulated  boxes)  were  swapped  from  one  site  to  another  on  several  occasions  as  service 
was  required.  To  draw  a  distinction  between  the  modem  sites  and  the  modem  packages, 
the  sites  were  numbered  1-6  and  the  packages  were  given  the  names:  bambi,  bashful, 
doc,  dopey,  gary,  wally.  The  data  collected  from  each  modem  is  stored  in  a  directory  of 
the  data  base  with  the  appropriate  package  name. 

The  location  of  the  modem  packages  throughout  the  9  data  capture  experiments  is  sum¬ 
marized  in  the  following  table.  A  blank  entry  indicates  that  no  modem  package  was 
present  at  the  site  concerned. 


TABLE  1.  Site  Table  (Experiment  Configuration) 


Experiment 

SKel 

Site  2 

sites 

Site  4 

Sites 

She  6 

l(rr) 

dopey 

bashful 

doc 

gary 

2  (mill) 

dopey 

bashful 

doc 

gary 

bambi 

3  (on27_3) 

dopey 

bashful 

gary 

bambi 

4  (on28_3) 

dopey 

bashful 

doc 

gary 

bambi 

5  (mit2) 

dopey 

bashful 

wally 

gary 

bambi 

doc 

6  (on28_3a) 

dopey 

bashful 

wally 

gary 

bambi 

doc 

7  (on29_3) 

dopey 

bashful 

wally 

gary 

bambi 

doc 

8  (on29_3a) 

dopey 

bashful 

wally 

gary 

bambi 

doc 

9  (collide) 

dopey 

bashful 

wally 

gary 

bambi 

doc 

4.3  Sources  of  Error  in  the  Data  Set 

A  number  of  error  and  interference  sources  are  known  to  have  influenced  the  experi¬ 
mental  data  set  in  terms  of  the  received  acoustic  signals,  the  reception  times,  the  GPS 
positions,  and  the  voltage  and  temperature  measurements.  These  are  summarized  in  the 
following. 

i)  Received  Acoustic  Signals:  Two  sources  of  noise/interference  were  observed  during 
the  acoustic  data  capture.  First,  the  ambient  noise  level  fluctuated  widely  depending 
upon  the  wind  speed.  At  high  wind  speeds,  considerable  impact  noise  from  blown  snow 
and  ice  particles  was  transmitted  from  the  ice  into  the  water.  This  noise  source  is  known 
to  have  a  substantial  component  in  the  frequency  range  used  for  acoustic  communica¬ 
tions  (lOkHz  -  20kHz).  An  hourly  log  of  Ae  wind  speeds  during  the  deployment  has 
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been  obtained  from  PMEL  at  NOAA  (Appendix  B)  and  an  attempt  will  be  made  to  cor¬ 
relate  the  wind  speed  with  the  ambient  noise  levels  in  the  captured  acoustic  data. 

The  second  interference  source  affecting  the  acoustic  data  was  the  acoustic  navigation 
beacons  used  by  the  MIT  group,  (jertain  of  these  beacohs  generated  Sms  pulses  with 
high  eneigy  in  the  communications  frequency  band.  The  pulses  occurred  every  second. 

ii)  Reception  Timing:  An  error  in  the  received  packet  timing  software  in  the  acoustic 
receiver  resulted  in  an  occasional  timing  ambiguity.  The  magnitude  of  the  ambiguity  is 
exactly  1  second  and  can  be  readily  corrected  in  the  data  without  loss  of  confidence.  It 
should  also  be  noted  that  if  the  timing  figure  of  merit  (TFOM)  of  the  GPS  is  >7  or  F,  the 
reception  timing  should  be  disregarded. 

Ui)  GPS:  The  accuracy  of  the  GPS  ‘xyz’  position  is  within  l(X)m  due  to  the  DOD  selec¬ 
tive  availability  randomization. 

iv)  Temperature  and  Voltage:  The  temperature  and  voltage  readings  logged  during 
Expts.  3, 4,6,7, 8  are  to  be  disregarded.  These  measurements  were  taken  with  an  essential 
power  supply  disabled  and  so  are  inaccurate. 
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5.0  Result  Formats 


5.1  Data  File  Format 

Data  collected  from  each  node  are  stored  in  files  which  are  named  with  the  following 
conventions;  out<minute>.<daldal>,  where:  out  is  the  prefix  for  all  output  data  files; 
<nunute>  is  the  minute  of  the  day  that  the  data  file  was  created  (computed  firom  GMT); 
and  <daldal>  is  the  file  extension.  For  example,  in  a  particular  experiment  subdirectory 
for  a  given  node  there  might  be  two  files  named:  6ut889.da  and  out889.dal.  Due  to  cer¬ 
tain  conditions  (copy  times  over  the  RF  link  and  security  circuit  time-out  period)  data 
had  to  be  transmitted  to  the  base  camp  as  two  separate  files.  Therefore,  a  file  with  a  dal 
extension  is  a  continuation  of  the  file  with  the  same  name  and  a  da  extension. 

Critical  information  about  the  usage  of  this  data  is  contained  in  the  header  of  the.da 
files.  This  information  has  the  following  format; 

1234  53352.(X)0000 

1234:  Arctic  Data  Magic  Number 
53352.0(XX)0()0:  Nearest  Second  of  Day  to  detect 
17120  17528 

17 120:  Index  into  data  for  start  of  packet 
17528:  Index  into  data  for  stop  of  packet 
43 

4;  Number  of  channels  of  data 

3:  Number  of  channels  on  which  a  detect  occurred 

There  is  one  entry  in  the  following  format  for  every  channel  on  which  data  was  col¬ 
lected. 

1  -0.693600 

1<0>;  1  is  detect  on  this  channel.  0  no  detect 

-0.693600:  Offset  in  seconds  from  Nearest  second  of  day  to  a  detect  to  the  actual 
detection. 

1  -0.693600 
1  -0.684200 
0-0.693600 

The  data  is  stored  in  two  columns,  the  left  hand  column  is  the  real  component  of  the 
data  and  the  right  hand  column  is  the  imaginary.  Data  samples  are  stored  in  ascending 
order  from  earliest  to  latest.  Multi-channel  data  is  stored  in  the  following  sequence: 

Channel  1,  Sample  1 
(Channel  2,  Sample  1 
Channel  3,  Sample  1 
Qiaimel  4.  Sample  1 
Channel  1.  Sample  2 
Channel  2.  Sample  2 
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5.2  STATUS.LOG  File  Format(s) 

There  are  two  STATUS.LOG  file  formats  as  indicated  below.  Both  contain  time  infor- 
matioi,  latitude,  longitude,  temperature  measurements,  battery  voltage  and  current  and 
an  indicator  as  to  whether  the  modem  was  transmitting  or  receiving  at  the  time. 

Format  A  is  used  in  experiments  1  through  4,  and  Format  B  is  used  in  experiments  6 
through  8.  Experiments  5  and  9  have  no  STAnJS.LOG  files. 

5.2.1  Format  A 

HH:MM:SS  TFOM  =  X  DDMM.HH.NJ)DDMM.HH 
TO  SDTO  T1  SDTl  T2  SDT2  T3  SDT3  BAIT  CURRENT 
txirxone 


HH;MM:SS  GMT  Hours  (HH):  Minutes  (MM):  Seconds  (SS) 

TFOM  =  X  Timing  Figure  of  Merit:  X  is  an  integer  between  1 

and  15  which  indicates  the  accuracy  of  the  timing 
pulse  from  the  GPS  (See  Table). 

DDMM.HH,NT)DDMM.HH  Latitude:  Degrees  (DD)  Minutes  (MM).  Hundredths 

of  Minutes  (HH),  North  (N),  Longitude:  Degrees 
(DDD)  Minutes  (MM).  Hundredths  of  Minutes  (HH) 
TO  SDT0...T3  SDT3*  (T0-T3)  Temperature  for  channels  0-3  in  degrees 

Celsius. 

(SDT0-SDT3)  Standard  Deviation  of  Temperature 
Measurement  for  channels  0-3  in  degrees  Celsius. 
BATT  Battery  voltage  (volts). 

CURRENT*  Supply  Cinrent  (amps). 

txirx  one  M^e  indicator:  (tx)  transmitter  mode,  (rx)  receiver 

one:  indicates  a  single  transmission  or  reception. 

*See  individual  experiment  description  for  details. 

5.2.2  R)rmat  B 

HH:MM:SSTFOM  =  X  DDMM.HHjNfJ)DDMM.HH 
TO  SDTO  T1  SDTl  T2  SDT2  T3  SDT3  BAIT  CURRENT 
XXXXYYYYZZZZ 
txirxone 


HH:MM:SS  GMT  Hours  (HH):  Minutes  (MM):  Seconds  (SS) 

TFOM  =  X  Timing  Figure  of  Merit:  X  is  an  integer  between  1 

and  15  which  indicates  the  accuracy  of  the  timing 
pulse  from  the  GPS  (See  Table). 

DDMM.HH,N  J)DDMM.HH  Latitude:  Degrees  (DD)  Minutes  (MM).  Hundredths 

of  Minutes  (HH),  North  (N),  Longimde:  Degrees 
(DDD)  Minutes  (MM).  Hundredths  of  Minutes  (HH) 
TO  SDT0...T3  SDT3*  (T0-T3)  Temperature  for  channels  0-3  in  degrees 
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BAIT 

CURRENT* 
txirx  one 

XXXXYYYYZZZZ 


Celsius. 

(SDT0-SDT3)  Standard  Deviation  of  Temperature 
Measurement  for  channels  0-3  in  degrees  Celsius. 
Battery  voltage  (volts). 

Supply  Current  (amps). 

Mode  indicator;  (tx)  transmitter  mode,  (rx)  receiver 
one:  indicates  a  single  transmission  or  reception. 

X,  Y  and  Z  coordinates  of  modem  in  meters. 


*See  individual  experiment  description  for  details. 


TABLE  2.  Timing  Figure  of  Merit  (TFOM) 


TFOM  Value 

Expected  Time 
Error  (min) 

Expected  Time 
Error  (max) 

1 

OnS 

InS 

2 

1-nS 

lOnS 

3 

lOnS 

lOOnS 

4 

lOOnS 

luS 

5 

luS 

lOuS 

6 

lOuS 

lOOuS 

7 

lOOuS 

ImS 

8 

ImS 

lOmS 

9 

lOmS 

lOOmS 

10-14 

NOT  USED 

NOT  USED 

15 

No  information 

Available 
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6.0  Experiment  Descriptions _ 

Details  of  the  9  data-logging  experiments  performed  are  given  in  the  following  sub-sec¬ 
tions.  ■ 

6.1  Experiment  1  (n) 

Description:  Round-robin  of  transmissions  from  each  modem  to  remaining 

modems. 

Data  purpose:  Analysis  of  network  connectivity,  modem  localization,  and 

comparative  performance  of  the  network  paths. 

Started:  8: 10,  March  26 

Duration:  approx.  3  hours 

Number  of  nodes:  4 

Number  of  channels:  1 

Packet  length:  1000  symbols 

Raw  data:  /dataOMarky/arctic/rr 

Processed  data:  /dataO/marky/arctic/exptl 


TABLE  3.  GPS  positions.  Expt  1 . 


Node 

GPS 

bashful 

72:56.50N.  148:26.95W 

doc 

72:57.04N,  148:29.91W 

dopey 

72:57.97N,  148:31.21W 

gary 

72:59.03N.  148:28.32W 
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6.2  Experiment  2  (mit1) 


Description: 

Data  purpose: 
Started: 

Duration: 

Number  of  nodes: 
Number  of  channels: 
Packet  length: 

Raw  data: 

Processed  data: 


Receptions  from  Datasonics  modem  on  the  Odyssey  (Odyssey 
approx.  8m  below  ice  surface). 

Analysis  of  AUV  up-link  performance,  AUV  localization. 
14:03,  March  27 
3  hrs  20  mins 
5 
1 

1000  symbols 

/dataOAnarky/arctic/mitl 

/data0/marky/arctic/expt2 


TABLE  4.  GPS  positions,  Expt  2. 


Node 

GPS 

bambi 

72:59.42N,  148:00.39W 

bashful 

72:56.86N.  147:58.87W 

doc 

72:57,39N,  148:0L83W 

dopey 

72:58.38N,  148:03.32W 

gary 

72:56.66N,  i48:13.33W 

Odyssey 

TBD 
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6.3  Experiment  3  (on27_3) 


Description: 
Data  purpose: 


Started: 

Duration: 

Number  of  nodes: 
Number  of  channels: 
Packet  length: 

Raw  data: 

Processed  data: 


Round-robin  with  hourly  change  of  transmitting  node. 

Analysis  of  long-term  netwOTk  performance,  acoustic  channel 
variability,  modem  localization,  reduced-complexity  receiver 
design. 

22:38,  March  27 
lOhrs,  23mins 
4 
1 

5000  symbols 

/data0/marky/arctic/on27_3 

/data0/marky/arctic/expt3 


TABLE  5. 


GPS  positions  and  experiment  times,  Expt.  3 


Node 

Start  Time 

GPS 

Stop  Time 

GPS 

bambi 

22:38:18 

7259.40 

14800.37 

08:39:59 

7258.91 

14755.66 

bashful 

22:38:20 

7256.85 

14758.84 

08:44:16 

7256.41^^1. 

14754.33 

dopey 

22:38:19 

7258.36JS[, 

14803.29 

08:47:02 

7257.87 

14758.57 

gary 

22:38:21 

7256.65 

14813.29 

08:43:23 

7256.23 JS[. 
14808.88 

TABLE  6.  Transmission  schedule,  Expt.  3 


Ttme 

Ihmsmitting 

Node 

File  Time  Stamps 

22:38-23:00 

bambi 

420-479 

23:00  -  0:00 

bashful 

480  -  539 

0:00  - 1:00 

gary 

540-599 

1:00-2:00 

no  transmitter 

600  -  659 

2:00-3:00 

dopey 

660-719 

3:00-4:00 

bambi 

720-779 

4:00-5:00 

bashful 

780-839 

5:00-6:00 

gary 

840  -  899 

6:00  -  7:00 

no  transmitter 

900  -  959 

7:00-8:00 

dopey 

960-1019 

8:00  -  8:40 

bambi 

1020-  1079 
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QPS  Tracking  of  (on27) 


FIGURE  1.  GPS  tracking,  Expt.  3. 
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6.4  Experiment  4  (on28_3) 


Description: 
Data  purpose: 


Started: 

Duration: 

Number  of  nodes: 
Number  of  channels: 
Packet  length: 

Raw  data: 

Processed  data: 


Round-robin  with  hourly  change  Of  transmitting  node. 

Analysis  of  long-term  network  performance,  acoustic  channel 
variability,  modem  localization,  reduced-complexity  receiver 
design. 

10:12,  March  28 
5  hrs,  10  mins 
3 
1 

10000  symbols 

/data0/marky/arctic/on28_3 

/dataOAnarky/arctic/expt4 


TABLE  7.  GPS  positions  and  experiment  times,  Expt.  4. 


Node 

Start  Time 

GPS 

Stop  Time 

GPS 

bashful 

10:12:20 

7256.40 JJ. 
14754.35 

15:21:58 

7256.39JSi, 

14754.49 

dopey 

10:12:13 

7257.91J^, 

14758.41 

15:20:33 

7257.87 

14758.67 

gary 

10:12:20 

7256.26 JT. 
14808.71 

15:20:05 

7256.20JN. 

14809.03 

TABLE  8.  Transmission  schedule,  Expt.  4 


Time 

lyansmitting 

Node 

File  Time  Stamps 

10:12-11:00 

dopey 

420  -  479 

11:00-12:00 

bashful 

480  -  539 

12:00  - 13:00 

gary 

540  -  599 

13:00  - 14:00 

dopey 

600  -  659 

14:00  - 15:00 

bashful 

660  -  719 

15:00  - 15:20 

gary 

720  -  779 
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6.5  Experiment  5  (mit2) 


Description: 

Data  purpose: 

Started: 

Duration: 

Number  of  nodes: 
Number  of  channels: 
Packet  length: 

Raw  data: 

Processed  data: 


Receptions  from  Datasonics  modem  on  the  Odyssey  (Odyssey 
.  approx.  8m  below  ice  surface). 

Analysis  of  AUV  up-link  performance,  AUV  localization, 
multi-channel  (diversity)  receiver  performance. 

22:00,  March  28 
20  mins 
6 
4. 

2000  symbols 

/data0Marky/arctic/mit2 

/dataOAnarky/arctic/exptS 


TABLE  9.  GPS  positions,  Expt  2. 


Node  GPS 


bambi 

72:58.78N,  147:57.08W 

bashful 

72:56.28N,  147:55.99W 

doc 

72:58.68N,  148;00.45W 

dopey 

72:57.74N,148:00.03W 

gary 

72:56.10N,  148:10.49W 

wally 

72:56.83N,  147:58.98W 

Odyssey 

TBD 
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6.6  Experiment  6  (on28_3a) 


Description: 
Data  purpose: 


Started: 

Duration: 

Number  of  nodes: 
Number  of  channels: 
Packet  length: 

Raw  data: 

Processed  data: 


Round-robin  with  hourly  change  of  transmitting  node,  multi¬ 
channel  reception. 

Analysis  of  long-term  network  performance,  diversity  recep¬ 
tion,  acoustic  channel  variability,  modem  localization, 
reduced-complexity  receiver  design. 

22:22,  March  28 
9  hrs,  50  mins 
6 
4 

2000  symbols 

/data0Anarky/arctic/on28_3a 

/data0/marky/arctic/expt6 


TABLE  10.  GPS  positions  and  experiment  times,  Expt.  6. 


Node 

Start  Time 

GPS 

Stop  Time 

GPS 

bambi 

22:20:27 

7258.88J^, 

14757.66 

08:27:47 

7259.73X 

14807.15 

bashful 

22:20:35 

7256.36J^, 

14756.48 

08:27:39 

7257.26X 

14806.16 

doc 

22:20:31 

7258.80 Jsr, 
14800.99 

06:40:39 

7259.33JSr, 
14808.20  • 

dopey 

22:20:17 

7257.843Sr, 

14800.58 

08:32:36 

7258.70JS[, 

14810.30 

gary 

22:20:25 

7256.18JSr, 

14811,03 

06:43:48 

7256.803Sf, 

14818.44 

wally 

22:20:22 

7256.91  JST. 
14759.50 

08:31:49 

7257.80J^, 

14809.43 

TABLE  11 .  Transmission  schedule,  Expt.  6 


Time 

Transmitting 

Node 

File  Time  Stamps 

22:20  -  23:00 

bambi 

420  -  479 

23:00  -  0:00 

dopey 

480  -  539 

0:00  - 1:00  • 

wally 

540  -  599 

1:00-2:00 

gary 

600  -  659 

2:00-3:00 

bashful 

660  -  719 

3:00-4:00 

doc 

720  -  779 

4:00-5:00 

bambi 

780  -  839 

5:00  -  6:00 

dopey 

840-899 

6:00-7:00 

wally 

900  -  959 

Report  on  the  Acoustic  Network  Arctic  Deployment,  March  1994 


19  of  50 


Experiment  Descriptions 


IVansmlttlng 

Time  Node  File  Time  Stamps 

7:00  -  8:00  gary  960-1019 

8:00  -  8:30  bashful  1020  - 1079 


QPS  Tracking  of  (on28a) 


HGURE  3.  GPS  tracking,  Expt.  6. 
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6.7  Experiment  7  (on29_3) 


Description: 
Data  purpose: 


Started: 

Duration: 

Number  of  nodes: 
Number  of  channels: 
Packet  length: 

Raw  data: 

Processed  data: 


Round-robin  with  hourly  change  of  transmitting  node,  multi¬ 
channel  reception. 

Analysis  of  long-term  network  performance,  diversity  equal¬ 
izer,  acoustic  channel  variability,  modem  localization,  reduced- 
complexity  receiver  design. 

9:57,  March  29 
11  hrs,  38  mins 
6 
4 

2000  symbols 

/data0Marky/arctic/on29_3 

/data0/marky/arctic/expt7 


TABLE  12. 

GPS  positions  and  experiment  times,  Expt.  7. 

Node 

Start  Time 

GPS 

Stop  Time 

GPS 

bambi 

09:57:39 

7300.00 JT, 
14809.26 

21:35:52 

7302.81  Jl, 

14825.08 

bashful 

09:57:51 

7257.54 Jl, 
14808.19 

21:30:07 

7300.30 JSf, 

14823.67 

doc 

09:57:35 

7259.90 JSr. 
14812.62 

21:06:30 

7302.38JSr, 

14826.99 

dopey 

09:57:48 

7258.94JJ. 

14812.20 

21:33:24 

7301.76J^, 

14828.07 

gary 

09:57:38 

7257.34 

14822.95 

21:13:45 

7300.15,N, 

14838.29 

wally 

09:57:53 

7258.093^, 

14811.39 

21:29:25 

7300.89 

14826.72 

TABLE  13. 

Transmission  schedule,  Expt.  7 

Time 

Transmitting 

Node  File  Time  Stamps 

9:57  - 10:00 

doc 

1080 

-1139 

10:00-11:00 

bambi 

1140 

-1199 

11:00-12:00 

dopey 

1200 

-1259 

12:00  - 13:00 

wally 

1260 

-1319 

13:00  - 14:00 

gary 

1320 

-1379 

14:00-15:00 

bashful 

1380 

-1439 

15:00  - 16:00 

doc 

0-59 

16:00  - 17:00 

bambi 

60-119 

17:00  - 18:00 

dopey 

120- 

179 
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lyansmittlng 

Node 


18:00  - 19:00  wally 
19:00-20:00  gary 
20:00  -  21:00  bashful 
21:00-21:29  doc 


File  Time  Stamps 

180  -  239 
240  -  299 
300-359 
360-419 


QPS  Tracking  of  (on29) 


O) 

Q  73- 


=  72.99- . 


Vn . 


****** ^ 


•\waly. 


bas 


gary_ 

72.95 1 - 1 - 1 - ^ - 1 - ^ - 1 

148.7  148.6  148.5  148.4  148.3  148.2  148, 

Longitude  (Degrees  E.) 


HGURE  4. 


GPS  tracking,  Expt.  7. 
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6.8  Experiment  8  (on29_3a) 


Description: 
Data  purpose: 


Started: 

Duration: 

Number  of  nodes: 
Number  of  channels: 
Packet  length: 

Raw  data: 

Processed  data: 


Round-robin  with  hourly  change  of  transmitting  node. 

Analysis  of  long-term  netwOTk  performance,  acoustic  channel 
variability,  modem  localization,  reduced-complexity  receiver 
design. 

22:40,  March  29 
9  hrs,  50  mins 
6 
1 

10000  symbols 

/data0/marky/arctic/on29_3a 

/dataO/marky/arctic/exptS 


TABLE  14.  GPS  positions  and  experiment  times,  Expt.  8. 


Node 

Start  Time 

GPS 

Stop  Time 

GPS 

bambi 

22:40:54 

7303.18X 

14826.41 

08:09:02 

7305.28 

14833.38 

bashful 

22:40:54 

7300.67 JS[. 
14825.16 

08:08:18 

7302.70J^. 

14832.04 

doc 

22:40:49 

7303.09 JSr. 
14829.77 

08:10:50 

7305.19X 

14836.73 

dopey 

22:40:52 

7302.14J^, 

14829.35 

08:11:02 

7304.24 JST, 
14836.33 

gary 

22:40:56 

7300.54 

14839.76 

08:10:35 

7302.56J1. 

14846.16 

wally 

22:40:51 

7301. 24 JM, 
14828.19 

08:07:36 

7303.29 J>I, 
14835.05 

TABLE  15.  Transmission  schedule,  Expt.  8 


Ttansmittlng 

Time  Node  File  Time  Stamps 


22:40  -  23:00 

bambi 

420  -  479 

23:00  -  0:00 

dopey 

480-539 

0:00-1:00 

wally 

540-599 

1:00-2:00 

no  transmitter 

600  -  659 

2:00-3:00 

bashful 

660  -  719 

3:00-4:00 

doc 

720  -  779 

4:00-5:00 

bambi 

780  -  839 

5:00-6:00 

no  transmitter 

840  -  899 

6:00  -  7:00 

wally 

900-959 
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Transmitting 

Time  Node  Fiie  Time  Stamps 

7:00-8:00  gary  960-1019 

8:00-8:10  bashful  1020-1080 


GPS  Tracking  of  (on) 


Longitude  (Degrees  E.) 


RGURE  5.  GPS  tracking,  Expt.  8. 
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6.9  Experiment  9  (collide) 


Description: 

Data  purpose: 
Started: 

Duration: 

Number  of  nodes: 
Number  of  channels: 
Packet  length: 

Raw  data: 

Processed  data: 


Near-simultaneous  transmissions  from  two  modems  (dopey 
and  bashful)  to  remaining  modems. 

Analysis  of  multi-user  and  full-duplex  receiver  algorithms. 
10:42,  March  30 
approx.  30  minutes 
6 
4 

2500  symbols 

/dataOAnarky/arctic/collide 

/data0/marky/arctic/expt9 


TABLE  1 6.  GPS  positions,  Expt  9. 


Node 

GPS 

bambi 

73:05.28N,  148:33.38W 

bashful 

73:02.70N.  148:32.04W 

doc 

73:05.19N,  148:36.73W 

dopey 

73:04.24N,  148:36.33W 

gary 

73:02.56N,  148:46.16W 

wally 

73:03.29N.  148:35.05W 
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7.0  Data  Extraction  and  Multi-channel  Receiver  Scripts 


7.1  MATLAB  script  for  automatic  gain  controi,  AGC.M 

function[y,p]  =  agc(x,  mu) 

% 

%[y,p]  =  agc(x,  mu) 

%Feedforward  automatic  gain  control  with  parameter  mu. 

%x  is  the  input  signal  vector  (time  series),  x  may 
%be  complex. 

%mu  sets  the  pole  position  of  the  first  order  input 
%power  filter  -  pole  is  at  1-mu,  A  small  mu 
%(e.g.,  0.005)  results  in  slow  age  action. 

%yis  the  output  signal  vector 
%pis  the  input  signal  power  estimate. 

% 


%Mark  JohnsonAcoustic  Telemetry  Group,  WHOI 
% 

%Last  Modified:24  October  1994 
% 


%  prime  power  estimate  with  a  power  estimate  from  the 
%  first  5  time  constants  of  x. 

nstart  =  ceil(5/mu) ; 

xs  =  abs(x).'^2 ; 

p_start  =  mean(xs(l:nstart)) ; 

start^gain  =  l/sqrt(p_start) ; 

%  compute  power  estimate  for  all  of  x 

p  =  filter([mu  0],  [1  -(1-mu)],  xs,  p^start) ; 

%  multiply  X  by  inverse  sqrt  of  power  estimate  to  get 
%  mean  power  of  1 

ageo  =  sqrt(p).^(-l) ; 
y  =  x.*agco ; 


7.2  MATLAB  script  for  demodulation  and  decimation,  DEMDEC.M 

function  Xdd  =  demdec(X,fcis,df) 

% 

%Xdd=demdec(X,fc,fs  ,df) 

%Demodulate  and  decimate  a  passband  signal  set. 

% 

%X  is  the  raw  input  data  sequence  in  column-signal 
%format. 

%fc  is  the  carrier  frequency. 

%fs  is  the  passband  sampling  frequency. 

%df  is  the  decimation  factor. 
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% 

%Xdd  is  the  demodulated  and  decimated  output  matrix  in 
%column-signal  format. 

%  . 

%  ' 

%  Mark  Johnson  Acoustic  Telemetry  Group,  WHOI 

% 

%  Last  Modified:  24  October  1994 

% 

[np  nch]  =  size(X) ; 

Xd  =  X.*(exp(-j*2*pi*fc/fs*(l:np))’)*ones(l,nch) ; 

Xdd  =  zeros(np,  nch) ; 
for  i=l:nch, 

Xdd(:,i)  =  decimate(Xd(:,i),df) ; 
end 


7.3  MATLAB  script  for  the  decision  feedback  equaiizer,  EQ.M 

function[W.de.dd.sef,F.fr]  =  eq(X,T4nap,pll,N,options,fcn) 

% 

%rW,de,dd,se,F,fr]  =eq(X,T,m  ap  ,plLf ,N,options,fcn) 

%Multichannel  joint  fractionally-spaced  RLS  and  PLL  equaliser  and 
%decoder  for  coherent  communications.  See  Stojanovic  et  al. 

%' IEEE  Trans. ... 

% 

%X  is  the  complex  input  data  sequence.  It  has  a  column 
%  for  each  observed  signal.  The  sampling-rate  for  all 
%channels  is  nnew  samples  per  symbol  (see  options  below). 

%The  algorithm  will  attempt  to  decode  Ns  symbols  where 
%Ns  =  (number  of  rows  in  X)  /  nnew. 

%T  is  the  complex  training  sequence.  The  initial  training 
%duration  is  determined  by  the  size  of  T 
%map  is  the  complex  alphabet  vector. 

%N  is  a  n+1  vector  where  n  is  the  number  of  columns  in  X. 

%The  first  n  elements  specify  the  number  of  parameters  in  the 
%feedforward  filter  associated  with  the  corresponding 
%column  of  X.  The  support  of  the  ith  feedforward  filter 
%is  thus  N(i)/nnew  symbols.  The  final  entry  in  N  specifies  the 
%number  of  parameters  in  the  feedback  filter.  As  the  decision 
%feedback  is  always  at  the  symbol-rate,  the  feedback 
%support  is  N(n+1)  symbols. 

%pil  sets  the  PLL  filter  coefficients:  the  first  row  of  pll_f  is 
%the  numerator  coefficients.  The  second  row  is  the  denominator. 

%  Coefficient  organisation  is  as  for  Matlab's  filter.m  function. 

%options  is  a  vector  of  options: 

%[nnew;  lambda;  dinit;  THRESH;  BUDGET;  Nrec] 

%where  (default  values  are  in  { }): 

%nnew  {2}  is  the  number  of  input  measurements  to  be  used 
%per  symbol.  1/nnew  is  the  fractional  spacing. 

%lambda  {0.995}  is  the  forgetting  factor  for  the  RLS  algorithm. 

%dinit  { 100 }is  the  initial  diagonal  value  of  the  inverse 
%  correlation  matrix  used  in  the  RLS  algorithm. 

%THRESH  {0}is  the  mse  level  in  dB  above  which  the  equaliser 
%  parameters  are  updated  using  the  RLS  algorithm 
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%after  initial  training. 

%BUDGET  {Ns} the  maximum  number  of  RLS  updates  permitted  in 
%decoding  X. 

%Nrec  {50} sets  the  period  in  symbols  between  parameter  value 

%reports  (see  discussion  of  W  below). 

%fcn  is  the  name  of  the  .m  file  to  be  used  for  equalizer  update 

%Available  update  scripts  are  'n2rls',  'nSrls',  'sqrrls'. 

% 

%W  is  a  matrix  each  column  of  which  is  the  RLS  parameter  vector 
%  recorded  every  Nrec  symbols. 

%de  is  the  Ns  vector  of  equaliser  outputs. 

%dd  is  the  Ns  vector  of  decision  ou^uts. 

%sef  is  the  Ns  vector  of  filtered  squared  errors. 

%F  is  the  Ns  X  n  matrix  of  the  input  phase  compensation  as 
%determined  by  the  phase-locked-loop  algorithm. 

%fris  the  fraction  of  decoded  symbols  for  which  an  RLS  update 
%was  performed,  fr  <-  BUDGET/Ns. 

% 


%  Mark  Johnson  Acoustic  Telemetry  Group,  WHOI 
%based  on  rls.m  (jointfsrls)  by  Milica  Stojanovic 
% 

%  Last  Modified:  31  October  1994 
% 


%  Initialisation... 

seLnum  =  [0.05  0] ;  %  squared-error  filter:  single  pole  @  0.95 

sef_denom  =  [1  -0.95] ; 

if  nargin<7 

fen  =  ’n2rls’ ;  %  default  RLS  update 

end 

if  nargin<6  %  set  default  options 

nnew  =  2 ; 
lambda  =  0,995  ; 
dinit  =  1 ;  ‘ 

THRESH  =  0; 

BUDGET  =  floor(size(X,l)/nnew) ; 

Nrec  =  50 ; 
else 

nnew  =  options(l) ; 
lambda  =  options(2) ; 
dinit  =  options(3) ; 

THRESH  =  options(4) ; 

BUDGET  =  options(5) ; 

Nrec  =  options(6) ; 
end 
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np  =  sum(N) ;  %  total  number  of  parameters 

Nt  =  max(size(T)) ;  %  number  of  training  symbols 

Ns  =  floor(size(X,l)/nnew) ;  %  total  number  of  symbols  to  decode 

no  =  size(X2) ;  %  number  of  feedforward  elements  (channels) 

IN  =  [0;cumsum(N)] ; 

Phi  =  2eros(np,l) ;  %  regressor  vector 

eval(['[Theta  Zeq]=’  fen  '(np,  dinit);’]) ;  %  initialize  RLS  algorithm 

update  =  [’[Theta  Zeq]='  fen  '(Phi,e,lambda,Theta2eq);'] ;  %  script  for  RLS  update 

nfill  =  max(N(l : no) ./nnew) ;  %  initially  hold-off  RLS  updates  for 

%  nfill  symbols  until  regressor  has 
%  some  data  in  it. 
j  =  sqrt(-l) ; 

[f  PLLs]  =  filter(pll(l.:),pll(2,:),0) ;  %  initialize  PLL  filter  states  and  output 

PLLs  =  PLLs*ones(l.no) ;  %  ...for  each  input  channel 

g  =  zeros(no,l) ; 

de  =  zeros(Ns,l) ;  %  record  of  equaliser  output 

sef  =  zeros(Ns,l) ;  %  record  of  filtered  squared  estimation  errors 

dd  =  zeros(Ns,l) ;  %  record  of  hard  decision  output 

F  =  zeros(Ns,no) ;  %  record  of  phase-locked  loop  adjustments 

W  =  zeros(np4oor(Ns/Nrec)) ;  %  record  of  equalizer  parameters 

[sef(l)  Zsef]  =  filter (sef_num,  sef_denom,  0) ;  %  initialize  squared-error  filter  states 

d  =  0  ;  %  0th  feedback  symbol  is  zero 

updateent  =  0  ;  %  zero  the  RLS  update  counter 

k=l; 

%  Decode  Ns  symbols... 
for  i=l:Ns, 

%  ...compute  data  symbol  estimate  and  estimation  error 

for  ii=  1  :no,  %  for  each  input  channel 

.  %  . .  .compute  phase-locked  loop  coefficient 

[f  PLLs(:,ii)]  =  filter(pll(l.:),pll(2,:).g(ii),PLLs(:,ii)) ; 

F(i,ii)  =  f ;  %  record  PLL  coefficient 

%  ...and  update  regressor 

Phi(IN(ii)+l:IN(ii+l))  =  [e3q)(-j*f)*X((i-l)*nnew+l:i*nnew,ii);  Phi(IN(ii)+l:IN(ii+l)-nnew)] ; 

%  ...apply  feedforward  filter 

p(ii)  =  Theta(IN(u)+l:IN(ii+l))’*Phi(IN(ii)+l:IN(n^  ; 
end 

Phi(IN(no+l)+l:np)  =  [d ;  Phi(IN(ncH-l)+l:np-l)] ;  %  add  last  symbol  as  feedback  to  regressor 
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t=  Theta(IN(no+l)+l:np)'*Phi(IN(ncH-l)+l:np) ; 
de(i)  =  sum(p)+t ; 

ifi>Nt, 

[ee  ind]  =  min(abs(map“de(i))) ; 

d  =  map(ind) ; 

else 

d  =  T(i); 
end 


%  apply  feedback  filter 
%  sum  partial  outputs  for  symbol  estimate 

.  %  determine  .ou4)ut  symbol 
%  if  not  training... 

%  ...find  nearest  symbol  to  de  in  map  vector 
%  if  training  pick  next  training  symbol 


dd(i)  =  d ;  %  record  ou^ut  symbol 

e  =  d  “  de(i) ;  %  compute  error 

%  filter  squared-error  and  record 

[sef(i)  Zsef]  =  filter(sef_num,  sef_denom,  e*e‘,  Zsef) ; 

%  check  if  RLS  update  is  required 

if((i>nfiU)  &  ((updatecnt  <=  BUDGET)  &  (sef(i)  >  THRESH))  I  (i<=Nt)) 

% . .  .compute  RLS  parameter  update 

eval(update) ;  %  do  update 

updatecnt  =  updatecnt+1 ;  %  ...and  increment  update  count 

end 

if(rem(i  J^rec)~0)  %  report  equalizer  parameters  every  Nrec  symbols 

W(:Jc)  =  Theta; 

j^rintfC%d  M) ; 

k  =  k+l  ; 

end 

%  ...compute  PLL  update 
g  =  imag(p.*conj(p-i-e)) ; 
end 

^rintf('Number  of  updates  =  %dVn'.  updatecnt) ; 
ff  =  updatecnt/Ns ; 

7.4  MATLAB  script  to  run  the  equalizer,  EQRUN.M 

%  eqrun.m 

%Example  script  for  running  eq.m  receiver  and  displaying  results 
%Assumes  input  variables  y  (baseband  signals)  and  ps  (code  detect  point) 

%obtained  using: 

%[y  ps]  =  prepinput(x,  code,  sym_rate) ; 

% 

%see  prepinput.m  for  further  details. 

% 

% 

%  Mark  Johnson  Acoustic  Telemetry  Group ,  WHOI 

% 
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%  Last  Modified:  31  October  1994 
% 

sta  =  ps  +  nnew*nullsym  ;%  start  point  in  signal  vector 

options  =  [nnew;  lambda;  dinit;  THRESH;  BUDGET;  50] ;  %  options  vector  for  eq.m 

nn  =  [N(l:size(y,2))  *nnew;  N(size(N,l))] ;  %  convert  feedforward  filter  symbol  supports 

%  as  given  in  rx_set  to  number  of  parameters 

[W  de  dd  sef  ph  fr]  =  eq(y(sta:sta+nnew*nsymbols-l,:),train(l:ntrain),map,pll,nn,options,eqfcn); 

figure(2); 

ne  =  plotresult(de,  dd,  sef,  ph,  train(l  :size(dd,l)),  ntrain) ;  %  plot  results  in  4-graph  format 


7.5  MATLAB  script  for  locating  barker  code  for  synchronization, 
FIND1STPEAK.M 

functionps  =  findlstpeak(dv,thresh) 

% 

%ps=findls4)eak(dv,thresh) 

%Find  the  first  peak  in  vector  dv  that  exceeds  a  threshold. 

% 

%dv  is  the  input  'decision'  sequence. 

%thresh  is  the  level  above  which  a  detection  is  made, 

% 

%ps  is  the  index  of  the  first  peak  in  dv. 

% 


%  Mark  Johnson  Acoustic  Telemetry  Group.  WHOI 

% 

%  Last  Modified:  25  October  1994 

% 

%  find  first  point  that  exceeds  the  threshold 
ps  =  min(find(dv>=thresh)) ; 

%  find  the  immediately  following  peak 

ps  =  ps  +  min(find(diff(dv(ps:size(dv,l)))<0))  - 1 ; 

7.6  MATLAB  script  for  RLS  updating  of  Equalizer  parameters,  N2RLS.M 

function[C  Jl]  =  n2rls(X,error,lambda,C.R) 

% 

%[C.R]  =  n2rls(X,eiTor,lambda,C,R) 

%Standard  n-squared  RLS  update  implementation 
%for  real  or  complex  signals. 

%X  is  the  regressor  vector  (nxl) 

%erroris  the  current  error  value  (i.e., 

%desired  -  CR) 

%lambdais  the  forgetting  factor  (0<lambda<l) 

%Choose  value  near  1  (e.g.,  0.995)  for 
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%slow  tracking.  Tracking  time  constant 
%is  approximately  1/(1 -lambda)  samples. 

%Cthe  current  parameter  vector  (nxl) 

%Rthe  current  inverse  auto-correlation  matrix 
%estimate  (nxn).  / 

% 

%[C,R]  =  n2rls(n,  dinit) 

%Initializes  C  with  zeros(n,l)  and  R  with  diag(dinit) 

% 

% 

%  Mark  Johnson  Acoustic  Telemetry  Group,  WHOI 

% 

%  Last  Modified:  24  October  1994 


if(nargin  ==  2) 
n  =  X; 
dinit  =  error ; 

R  =  dinit*eye(n) ; 

C  =  zeros(n,l) ; 
else 

Z  =  R*X; 

mu  =  l/real(lambda  +  X**Z) ; 
K  =  mu*Z; 

R  =  (R-K*Z’)/lambda; 

C  =  C-herror'*K; 
end 


7.7  MATLAB  script  for  plotting  the  input,  PLOTINPUT.M 

functionplotinput(x,y,pdre.ps,settings) 

% 

%plotinput(x,y,p,ire,ps,settings) 

%display  input  signal  properties  for  coherent  receiver. 

%To  be  used  with  prepinput.m.  See  prepinputm  for  details. 

%xis  the  complex  input  data  sequence.  It  has  a  column 
%  for  each  observed  signal. 

%yprocessed  signal 

%pestimated  input  signal  power  from  agc.m 
%psstart  point  in  x  found  using  findlstpeak.m 
%irechannel  impulse  response  estimate  from  cre.m 
%settings  settings  vector  (see  prqjinputm) 

% 

% 

%  Mark  Johnson  Acoustic  Telemetry  Group,  WHOI 
% 

%  Last  Modified:  31  October  1994 

% 


symrate  =  settings(i) ; 
intf  =  settings(2) ; 
ire_display  =  settings(6) ; 

clg 
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ns  =  size(x,l) ; 

ind  =  (0:ns-l)/(symrate*intf) ; 

subplot(221),  plot(ind,abs(x(:,l))),  grid,  title('raw  time  signal') 
xlabel(’seconds’),  ylabelCamplitude') 

subplot(222),  plot(ind,abs(y(:.l))),  grid,  titieCprocessed  time  signal') 
xlabel('secx)nds') ,  ylabel('amplitude') 

subplot(223),  plot(ind,p),  grid,  titie('Input  power  estimate') 
xlabel('seconds').  ylabel('power') 

nire  =  ceil(symrate*intf*ire__display) ; 

subplot(224i  plot(((0:nire-l)-ps)/(symrate*intf),  ire(l:nire)),  grid 

titleCChannel  impulse  response  estimate') 

xlabel('time  in  seconds') ,  ylabel('amplitude’)  * 

7.8  MATLAB  script  for  plotting  results,  PLOTRESULT.M 

function  ne  =  plotresult(de,dd,sef.ph.train,ntrain) 

% 

%ne  =  plotresult(de,dd,sef,ph,train,ntrain) 

%display  results  of  coherent  receiver,  eq,m 
%de  is  the  vector  of  equaliser  ou^^uts. 

%  dd  is  the  vector  of  decision  outputs. 

%  sef  is  the  vector  of  filtered  squared  errors. 

%  ph  is  the  matrix  of  the  input  phase  compensation  coefficients. 

%trainis  the  vector  of  transmitted  symbols 
%ntrainis  the  number  of  symbols  to  be  used  for  training 
% 

%neis  the  number  of  symbol  errors  in  the  decoded  packet 
% 

% 

%  Mark  Johnson  Acoustic  Telemetry  Group,  WHOI 
% 

%  Last  Modified:  31  October  1994 
% 


nsymbols  =  size(de,l) ; 
ne=n  ; 


clg 

subplot(221),  plot(10*logl0(sef)),  grid,  title('MSE') 
xlabelCtime  in  symbols'),  ylabel('dB') 

subplot(222).  plot(de(ntrain+l:nsymbols),'g+'),  grid 
axis([-2  2  -2  2]),  axis(axis),  title('Constellation') 
xlabel('real(de)'),  ylabel('imag(de)') 

subplot(223),  plot(ph),  grid,  title('Phase') 
xlabel('time  in  symbols'),  ylabelCradians') 

if  nargin  ==  6 

errors  =  abs(dd  -  train)  >=  1  ; 
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ne=  sum(errors(ntrain+l:nsymbols)) ; 

fprintf(’number  of  symbol  errors  after  training  =  %d\n’,  ne) ; 

subplot(224),  plot(errors),  grid,  titlef Error  distribution') 

xlabelCtime  in  symbols'),  ylabel('errors') 

end 

drawnow; 


7.9  MATLAB  script  for  preparing  the  input  file,  PREPINPUT.M 

function  [y,ps,p]  =  prepinput(x,code,settings) 

% 

%[y,ps,p]  =  prepinput(x,code,settings) 

%Prepare  input  signal  to  equalizer  by  (i)  hardlimiting, 

%(ii)  start  detection,  and  (iii)  age 
%xis  the  complex  input  data  sequence.  It  has  a  column 
%  for  each  observed  signal.  The  sampling-rate  for  all 
%channels  is  nnew  samples  per  symbol  (see  settings  below). 

%codethe  complex  start  code  at  a  sampling-rate  of  intf  (see 
%below)  input  samples  per  code  symbol, 

%settings  settings  vector  containing: 

%settings  =  [symrate  intf  limit  start_thresh  agc_mu  ire_display] 

%where: 

%symrateis  the  code  symbol-rate  in  Hz 

%intfis  the  number  of  input  samples  per  code  symbol 

%limitthe  hard  limit  threshold  as  a  percentage  of  input 

%events. 

%start_threshthreshold  above  which  a  start  may  be  detected. 

%See  findls^eak.m  for  details. 

%agc_mucontrols  pole  position  of  age  power  filter.  See 
%agc.m  for  details. 

%ire_displayhow  many  seconds  of  the  impule  response  estimate 
%to  display, 

% 

%yprocessed  signal 

%psstart  point  in  y  according  to  cre.m  and  findlstpeak.m.  ps  is  in 
%input  samples  not  symbols. 

%pestimated  input  signal  power  from  agc.m 
% 


%  Mark  Johnson  Acoustic  Telemetry  Group,  WHOI 
% 

%  Last  Modified:  31  October  1994 
% 


nsearch  =  1000 ; 

intf  =  settings(2) ; 
limit  =  settings(3) ; 
start_thresh  =  settings(4) ; 
agc_mu  =  settings(5) ; 

xl  =  rmoutliers(x,limit) ;%  hard  limit 

codedf  =  interp(code,intf) ;%  interpolate  start  code 
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ire  =  cre(xl(l:nsearch),  codedf,  1) ;%  cross-correlate  code  and  input 

[y  p]  =  agc(xl,  agc_mu) ;%  age  input  . 

dv  =  ire.4)(l:nsearch) ;%  decision  vector  =  cross-correlator 
%  output  /  input  power 

ps  =  findlstpeak(dv,  start_thresh)  ;%  find  start  peak  in  decicion  vector 
ifps~=  [] 

fprintf(’start  found  at  input  point  %dSn',  ps) ; 
else 

fprintf(’no  start  found  with  threshold  =  %dNn',  starLthresh) ; 
end 

figure(l) ;%  ...plot  results 
plotinput(x,  y,  p.  ire,  ps,  settings) ; 

7.10  MATLAB  script  for  converting  raw  data  to  compiex  vector,  RAW2VEC.M 

function  x=raw2vec(v) 

% 

%x==raw2vec(v) 

%Convert  2-column  format  raw  data  from  AMS  VI  .0  receiver 
%into  multi-channel  signal  matrix.  Supports  files  from 
%Lake  Winnepesaukee  (Feb.  1994)  and  Arctic  (March  1994) 

%expeditions. 

% 

%vis  the  raw  data  (a  2-column  real  matrix) 

%xis  the  complex  recovered  time  series  in  column- 
%channel  format 
% 

% 

%  Mark  Johnson  Acoustic  Telemetry  Group,  WHOI 

% 

%  Last  Modified:  24  October  1994 

% 

RXM_MAGIC=  1234; 

BUFSPACE=100; 
x  =  [] ; 

if(v(l,i)  ==  RXM_MAGIC) 

%  Arctic  vintage  data  with  circular  buffer  structure 
%  header  contains: 

% 

%magic  numbertime  (secs) 

%data  enddata  start 
%chl  detectchl  time  offset 
%  :  : 

%chn  detectehn  time  offset 

% 

%  remainder  of  file  is  [real  imag]  per  row 
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nch  =  v(3,l); 
pos  =  v(2.:)  +  1 ; 
n  =  size(v,l) ; 
st  =  nch  +  4; 

y  =  v(st:n,l)  +j*v(st;n,2) ; 
n  =  size(y,l) ; 

fprintf(’Arctic  data  file:  %d  channels,  %d  samplesNn',  nch,  n/nch) ; 

first  =  pos(2)-BUFSPACE*nch ; 
if(first<  1) 
first  =  n+first ; 
end 

if(pos(l)  =  1) 
y  =  y(first:n) ; 
else 

y  =  y([first:n  l:pos(l)-l]) ; 
end 

for  i=l:nch, 

x=  [x  y(i:nch:size(y.l))] ; 
end  . 

else 

%  raw  data  is  Winn^esaukee  vintage 

fprintf(’Pre- Arctic  single-channel  data  fileNn') ; 
pos  =  v(i,:)+l ; 
n  =  size(v,l) ; 

X  =  v(2:n,l)  +  j*v(2:n,2) ; 
if(pos(l)  =  l) 
x  =  x(pos(l):n-l) ; 
else 

x  =  x([pos(l):n-l  l:pos(l)-l]) ; 

end 

end 

7.11  MATLAB  script  for  removing  outliers  from  data  set,  RMOUTLIERS.I\/l 

function  y  =  rmoutliers(x,  thresh) 

% 

%y  ss  rmoutlLers(x,  thresh) 

%Remove  outliers  from  time  series,  x,  using  symmetric 
%hard  limiter  with  limit  thresholds  set  from  the 
%histogram  of  x. 

% 

%x  is  the  signal  in  column  vector  format  (maybe 
%real  or  complex) 

%thresh  is  the  clipping  threshold  in  percentage  of 
%input  signal  events,  e.g.,  thresh  =  0.01  sets 
%the  limiting  level  such  that  only  0,01 

%of 
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%input  points  will  be  clipped. 

% 

%y  is  the  output  column  vector. 

% 

%  ■ 

%  Mark  Johnson  Acoustic  Telemetry  Group,  WHOI 

% 

%  Last  Modified:  24  October  1994 

% 

nbins  =  50  ;%  resolution  of  histogram 
ax  =  abs(x) ; 

[n,  xdist]  =  hist(ax,  nbins); 

lt  =  xdist(max(find(100*(l-cumsum(n)/sum(n))  >  thresh))) ; 
y  =  (ax  <=  lt).*x  +  lt*(ax  >  lt).*sign(x) ; 

7.12  MATLAB  script  with  example  settings  for  single  channel  reception, 
TEST_SET.M 

%rx_set.m 

%standard  receiver  settings  for  eq.m  and  prepinputm 
% 

% 

%  Mark  Johnson  Acoustic  Telemetry  Group ,  WHOI 

% 

%  Last  Modified:  31  October  1994 
% 


load  train  1  %  training  symbols 

%  settings  for  prepinput.m 

symrate  =  2500  ;  %  symbol  rate  (Hz) 

limit  =  0,04;  %  clipping  threshold  in  %  of  events 

agc^mu  =  0.005  ;  %  sets  pole  position  of  age  power  filter 

%  -  pole  is  at  l-agc_mu 

starLthresh  =  100 ;  %  power  ratio  threshold  for  detector 

ire_display  =  0.02  ;  %  maximum  length  of  impulse  response  to  display  (secs) 

intf  =  2  ;  %  number  of  samples  per  code  symbol  in  input 

%  ...settings  vector  for  prepinput.m 

code  =  M;l;-l;l;-l;-l;l;l;-l;-l;-l;-l;-l] ;  %  13-element  barker  start  code 

settings  =  [symrate  intf  limit  start_thresh  agc_mu  irc_display] ; 


%  settings  for  eqrun.m 

N  =  [5;9] ;  %  equalizer  feedforward  and  feedback  support 

%  in  symbols. 

map  =  [l+j;-l-j;-l+j;l-j];  %  map  vector  for  QPSK 

pll  =  [0.011  -0.01  0  ;  1  -2  1]  ;  %  pll  filter  numerator  and  denominator 
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%  ...this  pll  has  two  integrators  and  one  zero. 

nuUsym  =  20 ;  %  number  of  symbols  between  start  code  and  data 

THRESH  =  10^(-7.5/10) ;  %  mse  threshold  in  dB  -  equalizer  parameters  are  only 

%  adjusted  when  the  current  mse  is  above  THRESH 

nsymbols  =  1000  ;  %  number  of  symbols  to  decode 

nnew  =  2 ;  %  number  of  samples  per  symbol  in  input 

BUDGET  =  nsymbols  ;  %  maximum  number  of  RLS  updates  to  perform 

ntrain  =  200  ;  %  number  of  training  symbols  at  start  of  packet 

lambda  =  0.995  ;  %  RLS  foigetting  factor 

dinit  =  1000  ;  %  RLS  initial  regressor  power  value 

eqfcn  =  ’n2rls' ;  %  RLS  update  function  to  use. 

7.13  Source  Code  for  extracting:  Time,  Latitude  and  Longitude  from 


STATUS.LOG  files: 

- */ 

/*STATFILTC  -  This  code  extracts  time,  latitude  and  longitude  from  */ 

/*  STATUS.LOG  files  created  during  the  Arctic  deployment  of  March  1994  */ 

/*  The  program  removes  any  measurements  which  were  taken  with  a  TFOM  >10  */ 

/*  See  Magellan  data  sheet  for  details  oaTFOM  */ 

/*  */ 
f*  Author:  D.  Herold,  WHOI  V 

f*  Date:  4/26/94  */ 

/* - */ 


#include  <stdio.h> 

#include  <ctype.h> 

#include  <stdlib.h> 

#define  CORRECT^^READS  3 

main(aigc,argv) 

int  argc; 
char  *argv[]; 

{ 

FILE  *fnamel; 

FILE*fname2; 
char  in_str[80]; 
inti; 

inttfom; 

inthour,minute,second; 
floatiatitude,longitude; 

if  (  argc  !=  3 )  /*  check  for  correct  usage*/ 

{ 

fprintf(stdout,"Usage:  statfilt  <filel>  <file2>\n"); 
fprintf(stdout,"NOTE:  <filel>  must  be  a  dos2unix'd  status.log  file....\n"); 
exit(O); 

} 

/*  open  files*/ 


/*  file  pointers*/ 

/*  input  string*/ 
/*  declare  vars.*/ 
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if  ( (fnamel  =  fopen(argv[l],"r"))  ==  NULL ) 

{ 

exit(O); 

} 

if  ( (fname2  =  fopen(argv[2]  ,"w"))  ==  NULL ) 

{ 

exit(O); 

} 

/*  read  file  and  get  pertinent  information,  store  in  file2*/ 

while  ( fscanf(fnamel,”%s",in_str)  !=  EOF ) 

{ 

if  (CORRECT_READS  ==  sscanf(in_str,"%2d:%2d:%2d  ",&hour,&minute,&second) ) 

{ 

for  ( i  =  0;  i  <=  2;  i++) 
fscanf(fiiamel  ,”%s",in_str); 
sscanf(in_str,”%x"  ,&tfom) ; 
if  (tfom  <  10) 

{ 

fscanf(fiiamel,"%s”in_str): 
sscanf(in_str,"  %f ,N,%f '  ,&latitude.&longitude) ; 
fprintf(fname2,"%2d  %2d  %2d  ",hour,minute,second); 
l^rintf(f name2  ."%7 . 2f  %8 . 2fsn''  .latitude, longitude) ; 

} 

} 


} 

fclose(fnamel);/*  close  files*/ 
fclose(fiiame2); 

}  /*  exit*/ 


7,14  MATLAB  Script  for  for  correcting  and  plotting  GPS  Lattitude  and 
Longitudes  extracted  from  STATUS.LOG  files. 

%  gpsplot.m  D.  Herold,  WHOI  26  April  1994 
%  This  file  loads  data  from  an  Arctic  statusJog  file  which 
%  has  been  dos2unixed  and  then  run  though  the  statfilt.exe 
%  program  in  ^/aalan/arctic/gpsdata/  directory. 

% 

%  function  [pos.samples.etime]  =  gpsplot(experiment) 

% 

function  [pos.samples.etime]  =  gpsplot(experiment) 

nodes  =  [’bambj 
'bashj 

'doc _ ’ 

'dopyj 

'garyj 
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’walyj] ; 

nodenum  =  0 ; 
x_min  =  999 ;  . 
x_max'= -999 ; 
y_min  =  999 ; 
y_max  =  -999 ; 

pos  =  zeros(150,12) ; 
etime  =  zeros(150»6) ; 

clg 

for  k  =  l:length(nodes) 

file  =  [nodes(k,:)  experiment  '_gps’] ; 
filename  =  [file  Mogl ; 
if  exist(filename) 
filename 

nodenum  =  nodenum  +  1 ; 
load(filename) ; 
eval([’y=  file’;’]); 

%  get  the  elapsed  e?q)eriment  time 

tm  =  y(:,i).*60+y(:,2)+y(:,3)./60 ; 

tm(find(tm=max(tm))+l:length(tm))  =  tm(find(tm~max(tm))+l:length(tm)) 
+  24*60; 

etime(l:length(tm)4iodenum)  =  tm ; 

%  Convert  to  decimal  degrees 

y(:,4)  =  fix(y(:,4)./100)+(rem(y(:,4)./100,l)./60)  *100 ; 

pos(l:length(y).2*nodenum-l)  =  y(:,4) ; 

y(:,5)  =  fix(y(:.5)7100)+(rem(y(:,5)./100,l)./60)  *100  ; 

pos(l:length(y),2*nodenum)  =  y(:,5) ; 

samples(nodenum)  =  length(y) ; 

index(nodenum)  =  k; 

%  Get  max  and  min  for  axes 

x_min  =  min(x_min,min(pos(l:samples(nodenum).2*nodenum))) ; 
x_max  =  max(x_max,max(pos(l:samples(nodenum),2*nodenum))) ; 
y_min  =  min(y_min,min(pos(l:samples(nodenum),2*nodenum  - 1))) ; 
y_max  =  max(y_max,max(pos(l:samples(nodenum),2*nodenum  - 1))) ; 
else 

k  =  k+l  ; 
end 

end 

x„min  =  x_min  -  0.01 ; 
x_max  =  x_max  +  0.01 ; 
y„min  =  y_min  -  0.01 ; 
y_max  =  y_max  +  0.01 ; 
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delta^lat  =  y_max-y_min ; 

x_max  =  x_inin  +  deltajat  /  cos(y_inax*(pi/180)) ; 

figure(l) ;  %  plot  out  gps  positions  of  nodes 

axis('square'); 

axis('equal’); 

axis([x_min  x_max  y_min  y__max]); 
fork=  linodenum 

plot(pos(l:samples(k),2*k),pos(l:samples(k),2*k-l)/.’),  grid  on 
hold  on 

plot(pos(  1 ,2*k)  ,pos(  1 ,2*k- 1)  /x') 

text(pos(1.2*k)*0.02,pos(l,2*k-l)-0.002,nodes(index(k),:)) 

end 

view( 180,-90) 

Title  =  ['GPS  Tracking  of  (’  ejqjcriment ')'] ; 
title(Title) 

xlabelCLongitude  (Degrees  E.)') 
ylabelCLatitude  (Degrees  N.)’) 
hold  off 

figure(2)  %  not  required,  however,  lets  look  at  times  of  experiments 
for  k  =  Imodenum 

plot(etime(l:samples(k)  Jc)-max(min(etime)),’.'),  grid  on 
hold  on 

text(samples(k)  ,etime(samples(k)  Jc)-max(min(etime))  ,nodes(index(k), :)) 
end 

Title  =  ['GPS  Time  vs.  Sample  of  ('  experiment ')'] ; 

title('Iitle) 

xlabel('Sample') 

ylabel('Sample  Time  (Minutes)') 
hold  off 

figure(l) 
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Appendix  A  -  Log  Book  Notes _ 

The  following  information  is  extracted  from  the  log  book  kept  during  the  experiment. 

0.1  March  17, 1994 

•  1700  MJ,  DH  Arrive  Deadhorse,  AK. 

•  1900  Move  equipment  into  lab  space. 

0.2  March  18, 1994 

•0730  Move  all  equipment  into  new  lab  space  and  unpack. 

• 0900  Hook  up  TELEM  (SPARC  10),  restore  new  data  from  backups  of  work  done 
at  WHOI  after  air  shipment.  Having  the  same  problem  with  the  DOS  shell 
program.  What  you  type  in  the  dos  window  is  not  displayed  on  the  screen. 
Called  BE  at  WHOI.  There  was  no  easy  way  to  fix  the  problem.  Decided  to 
use  JOSKO  (SPARC  IPO  as  we  had  in  Winnipesaukee  experiment. 

•1130  Lunch 

•  1200  Finished  switching  IPC  in  for  SPARC  10 

•  1500  Sneezy,  Wally,  Doc  all  Tx  and  Rx  ok.  Hook  up  GPS  to  Sneezy  to  see  if  we  get 

reception  inside  builing. 

•  1615  Start  bringing  in  PVC  tubes  for  thru-ice  fittings  to  warm  them  up  prior  to 

assembly. 

0.3  March  19, 1994 

•  0800  No  journal  entries  for  this  day.  Presume  we  spent  it  chasing  apparent  software 

bug. 

0.4  March  20, 1994 

•  17(X)  JC  arrives  Deadhorse 

•  2000  Determine  that  problems  which  have  looked  like  transmitter  and  receiver  soft¬ 

ware  problems  are  actually  hardware  problems  not  previously  encountered  on 
A/D  board.  We  have  to  date  been  unlucky  and  only  tried  combinations  of 
modems  with  bad  boards. 

0.5  March  21, 1994 

• 0200  Finish  debugging  and  preparing  six  of  the  modems  for  deployment  initially  on 
ice.  Pack  equipment  for  trip  to  ice. 

•  0900  Depart  Deadhorse  for  Ice  Camp  East. 

•  1100  Arrive  Ice  camp,  unload  plane,  move  equipment  to  our  ice  hut. 

•  13(X)  Short  lunch  break 
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•  1400  Raise  base-camp  antenna,  organize  hut,  get  base  camp  electronics  running.  We 

are  having  problems  with  the  UPS,  we  keep  getting  “power  hits”.  We  assume 
that  it  is  due  to  us  overloading  the  UPS. 

0.6  March  22, 1994 

•  0800  It  turns  out  that  the  camp  generator  does  not  have  a  closely  regulated  frequency 

and  it  is  too  difficult  for  the  UPS  operate  correctly. 

•  0830  Deployed  Doc  -'SOm  from  camp:  GPS  position  7257. 16,N,  14822.90  Must  pre¬ 

pare  pvc  tube  endcaps  inside  after  they  have  wanned  up  inside  the  huts.  Do 
NOT  use  teflon  pipe  sealer. 

•  1300  Deploy  Bambi,  GPS  7258.04,N,  14827.70 

0.7  March  23, 1994 

•  1300  Deploy  Bashful,  GPS  7256.58,N,14823.83 

0.8  March  24, 1994 

0.9  March  25, 1994 

•0900  DH  departs  Ice  camp.  ED  arrives. 

0.10  March  26, 1994 

0.11  March  27, 1994 

•  0930  Did  temperature  log.  Bambi  taken  by  helicopter  to  replace  gary  who  wend 

down  yesterday,  onite  batch  file  running  on  dopey,  gary  and  doc,  failed  after 
first  iteration.  Made  some  changes  to  hrt.bat  (ovemite  batch  file)  adn  renamed 
collect.bat.  This  needs  the  following  environment  variables:  TXSEQ  (>=0.,  < 
NNODES).  NNODES  e.g.  4.  PRESLEEP  (sleep  time  between  Rx  and  copy). 
RXSLEEP  (sleep  time  in  seconds  between  Rx’s).  TXSLEEP  (sleep  time  in  sec¬ 
onds  beteen  Tx’s.  DATADIR  (destination  data  directory  on  e:.  These  are  all 
set  in  onset.bat  which  is  called  from  autorun.bat.  Also  m  mailbox.bat  @  start 
there  is  a  checd  for  e:%sysdir/collect.mbx,  only  if  it  exists  can  you  call  collect. 

•  1000  Gary  replaced  by  Bambi.  Bambi  came  up  good.  New  GPS  fix: 

7259.65,Na4814.10.  Temperatures:  17.4,25.5-14.3,-22.7 

•  1030  File  set  required  for  collect  is:  collect.bat,  onsetbat,  autorun.bat,  mailbox.bat, 

lsleep.bat,  remain.bat,  rx_mess.bat 

•  1030  Qiecked  LSLEEP  with  env.  variables,  works  ok. 

•  1215  Foung  problem  with  gary  was  (suspected)  faulty  ARLAN.  Repalced  with 

ARLAN  from  WaUy.  Also,  deadperson  battery  was  flat,  replaced  and  checked 
thru  checklist  ok.  Took  gary  to  new  site  with  helicopter.  Able  to  ARLAN 
communicate  ok.  GPS  is  7257.1 0,N,  14824.42 
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•  1403  MIT  1  data  set  taken? 

•  1723  Mm  data  set  done? 
•2236  Temperature  log 
•2240  Run  COLLECT 


0.12  March  28, 1994 

•0845  Stopped  COLLECT  -  all  except  Doc  returned.  Also  doc  did  not  collect  any 
paketes.  All  data  and  logs  are  in :  %datadir%/on27_3. 


TABLE  1 7.  Transmitters  for  2240  COLLECT 


GPS  Hour 

lYansmitter 

Time 

Stamp 

7 

bambi 

420479 

8 

bashful 

480-539 

9 

gary 

540-599 

10 

600-659 

11 

dopey 

660-719 

12 

bambi 

720-779 

13 

bashful 

780-839 

14 

gary 

840-899 

15 

900-959 

16 

dopey 

960-1019 

17 

bambi 

1020-1079 

•  1446  Datasonics  modem  on  Odyssey  is  doing  occasional  broadcasts  for  next  two 

hours  -  may  be  received  by  the  long  packet  capture  on  gary,  dopey,  and  bashful 
(on28_3) 

•  1500  Field  service  on  bambi  (copied  new  auitorun.bat)  installed  doc  in  new  site. 

•  15(X)  histalled  wally  in  camp  site  (was  doc  site) 

•1610  GPS  log 

•  1700  Temperature  log 
•2100  Temperature  log 
•2130  GPS  log 
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•2200  GPSgeLxyz. 


TABLE  18.  XYZ  Positions 


NODE 

X 

Y 

Z 

doc 

-1588321.42 

-991950.08 

6076456.88 

bambi 

-1587216.75 

-993415.57 

6076521.02 

dopey 

-1589628,87 

-993049.73 

6075957.36 

wally 

-1590715.19 

-994423.98 

6075441.62 

gary 

•  -1595155.84 

-989769.79 

6075038.56 

bashful 

-1590679.12 

-996339.06 

6075141.54 

•2200  captured  2  data  sets  from  modem  on  odessey.  These  are  4  chaimel  IK  sym¬ 
bols  and  are  stored  in  Mrr2  directories. 

•  2222  COLLECT  started  with  4  channel  rxdl .m 

0.13  March  29, 1994 

•  0830  Stopped  COLLECT  (doc  and  gary  had  reset  earlier  '^1540  rather  than  1730). 

•  0830  Temperature  log: 

•  0830  GPS  positions: 

•  0830  GPS  geCxyz: 

•0930  Re-start  COLLECT  with  multi-channel  short  packets.  Last  nights  collect 
results  are  stored  in  ON28_3_A  directories.  Transmit  schedule  was: 


TABLE  19. 


Transmitter  schedule 


GPS  Hour 

7 

8 

9 

10  • 

11 

12 

13 

14 

15 

16 
17 


transmitter 

bambi 

dopey 

wally 

gary 

bashful 

doc 

bambi 

dopey 

wally 

gary 

bashful 


OUT  Keys 

420:479 

480:539 

540:559 

600:659 

660:719 

720:779 

780:839 

840:899 

900:959 

960:1019 

1020:1079 
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•  1705  Stopped  data  capture.  AH  units  returned. 

•  1705  Temperature  log. 

•1711  GPS  Positions. 

•  17 17  Re-started  COLLECT,  (short  packet,  4  channel)  with  reduced  PRESLEEP, 

RXSLEEP  increments  (20  of  25  previously)  and  reduced  TXSLEEP  (40  of  50 
previously). 

•  1845  MIT  doing  Odyssey  ops  with  datasonics  modem  on  (1  per  minute). 

•  1915  Mir  stops  Odyssey. 

•2130  COLLECT  stopped  (multichaimel,  short  packets)  stored  in  ON29_3  directory. 
•2130  GPS  log 
•2130  Temperature  1(^ 

•  2245  COLLECT  started  with  single  channel  long  packets.  Previous  collect  stored  in 

ON29_3  directory. 

0.14  March  30, 1994 

•  1005  Flow  is  cracking.  ED  recovering  nodes  electronics.  Last  nites  collect  is  in  ON 

directories  (status.log’s  are  still  <mi  PCMCIA)  Did  some  collection  for  milica: 
dopey  and  bashful  transmitting,  remainder  receiving.  Data  is  4  channel,  3000 
symbols  (2500  transmitted),  stored  in  root  directories.  e.g.  gary_d. 

•  1015  TARd  directwies  to  tapes. 
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Appendix  B  >  Wind  Speeds 


The  wind  measurement  station  was  by  the  helo  hut  until  about  5  or  6  on  29  March,  so 
the  early  wind  directions  aren’t  too  good  but  should  be  in  the  right  ballpark.  Tunes  are 
GMT.  The  station  should  have  taken  one  sample  per  hour,  but  its  timer  ran  roughly  a 
minute/hour  slow(sometimes  worse).  The  10-minute  sample  interval  was  controlled  by 
a  more  accurate  clock,  so  that  windspeed  calculations  should  be  OK.  \\Tnd  direction 
is  TO,  not  from,  kind  of  trae  (magnetic+30  degrees).  Floe  rotation  is  already  included 
in  wind  direction. 


YrMoDaHrMi 

AirTCO 

WSpd(m/s)  WDir 

CompCFloe  Rotation) 

9403270515 

-27.36 

1.76 

280.00 

22.40 

9403270618 

-26.40 

3.04 

274.40 

21.00 

9403270720 

-25.86 

2.40 

271.60 

21.00 

9403270821 

-23.34 

0.80 

273.00 

21.00 

9403270923 

-25.38 

2.40 

259.00 

21.00 

9403271025 

-27.96 

0.96 

235.20 

22.40 

9403271127 

-28.20 

2.88 

27.80 

22.40 

9403271231 

-29.76 

1.12 

50.20 

22.40 

9403271331 

-30.54 

1.44 

257.60 

22.40 

9403271432 

-30.06 

2.24 

329.00 

22.40 

9403271535 

-30.42 

2.88 

267.40 

22.40 

9403271739 

-29.22 

7.04 

18.00 

19.60 

9403271841 

-28.14 

6.88 

260.40 

21.00 

9403271943 

-27.42 

7.84 

257.60 

19.60 

9403272044 

-26.82 

7.20 

266.00 

18.20 

9403272148 

-26.52 

7.52 

266.00 

18.20 

9403272248 

-26.34 

8.96 

266.00 

19.60 

9403272350 

-26.52 

8.32 

270.20 

18.20 

9403280054 

-26.94 

11.52 

268.80 

18.20 

9403280154 

-27.12 

6.56 

270.20 

18.20 

9403280256 

-27.96 

3.52 

274.40 

19.60 

9403280358 

-28.98 

3.68 

274.40 

19.60 

9403280459 

-29.58 

5.28 

273.00 

19.60 

9403280502 

-29.58 

5.28 

273.00 

19.60 

9403280602 

-30.00 

528 

275.80 

21.00 

9403280704 

-30.30 

2.72 

282.80 

21.00 

9403280809 

-30.60 

2.24 

278.60 

21.00 

9403280909 

-30.84 

2.24 

281.40 

21.00 

9403281009 

-31.14 

2.24 

292.60 

21.00 

9403281113 

-31.38 

2.40 

302.40 

21.00 

9403281213 

-31.74 

1.60 

281.40 

21.00 

9403281315 

-31.98 

1.44 

294.00 

21.00 

9403281417 

-32.40 

1.60 

277.20 

19.60 

9403281519 

-32.04 

1.12 

278.60 

21.00 

9403281621 

-31.20 

0.96 

277.20 

18.20 

9403281723 

-31.20 

0.96 

277.20 

18.20 

9403281825 

-30.30 

1.12 

303.80 

19.60 
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YrMoDaHrMi  AirT(C)  WSpd(m/s)  WDir  C<»np(Roe  Rotation) 


9403281926 

-29.52 

0.48 

284.20 

18.20 

9403282028 

-28.20 

0.96 

305.20 

16.80 

9403282130 

-28.26 

2.56 

259.00 

16.80 

9403282232 

-27.96 

2.56 

260.40 

16.80 

9403282334 

-27.30 

1.92 

256.20 

16.80 

9403290047 

-20.52 

15.68 

259.00 

16.80 

9403290552 

-26.40 

2.88 

287.00 

58.80 

9403290654 

-25.74 

3.52 

294.00 

58.80 

9403290755 

-24.96 

4.32 

294.00 

58.80 

9403290859 

-24.18 

4.64 

287.00 

58.80 

9403290959 

-23.10 

5.12 

282.80 

58.80 

9403291000 

-23.10 

5.12 

282.80 

58.80 

9403291101 

-22.14 

6.56 

289.80 

57.40 

9403291203 

-21.42 

8.48 

289.80 

57.40 

9403291305 

-21.60 

8.96 

292.60 

57.40 

9403291409 

-22.20 

10.72 

289.80 

57.40 

9403291509 

-21.42 

10.56 

288.40 

57.40 

9403291610 

-20.88 

17.44 

289.80 

57.40 

9403291713 

-20.10 

16.16 

292.60 

57.40 

9403291814 

-19.20 

15.52 

292.60 

57.40 

9403291917 

-18;48 

15.04 

288.40 

57.40 

9403292017 

-17.88 

15.84 

292.60 

57.40 

9403292120 

-17.34 

18.72 

296.80 

57.40 

9403292221 

-16.56 

20.96 

296.80 

57.40 

9403292323 

-15.96 

19.04 

294.00 

57.40 

9403300027 

-15.78 

21.76 

299.60 

57.40 

9403300127 

-15.30 

18.56 

296.80 

57.40 

9403300228 

-15.06 

15.36 

302.40 

57.40 

9403300331 

-14.82 

17.28 

302.40 

57.40 

9403300431 

-14.70 

20.16 

305.20 

57.40 

9403300534 

-15.84 

16.80 

309.40 

58.80 

9403300635 

-15.84 

16.80 

309.40 

58.80 

9403300738 

-16.56 

10.40 

309.40 

58.80 

9403300841 

-17.16 

11.68 

305.20 

58.80 

9403300941 

-17.40 

11.84 

309.40 

58.80 

9403301045 

-17.52 

9.60 

309.40 

58.80 

9403301145 

-17.70 

14.24 

312.20 

58.80 

9403301246 

-17.94 

7.84 

312.20 

58.80 

9403301349 

-17.82 

6.56 

310.80 

58.80 

9403301449 

-17.64 

7.04 

317.80 

58.80 

9403301554 

-17.16 

7.04 

319.20 

58.80 

9403301654 

-17.04 

5.28 

324.80 

58.80 

9403301755 

-16.50 

4.32 

330.40 

57.40 

9403301859 

-16.20 

3.52 

350.00 

57.40 

9403301959 

-15.42 

2.40 

6.80 

58.80 

9403302001 

-15.42 

2.40 

6.80 

58.80 

9403302101 

-15.00 

1.60 

313.60 

57.40 
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YrMoDaHrMi  AirT(C)  WSpd(m/s)  WDir  Comp(Floe  Rotation) 


9403302203 

-14.88 

2.08 

340.20 

58.80 

9403302304 

-14.22 

2.72 

62.80 

58.80 

9403310008 

-14.58 

2.88 

58.60 

58.80 

9403310108 

-14.46 

3.20 

48.80 

58.80 

9403310209 

-14.40 

4.16 

296.80 

58.80 

9403310313 

-14.34 

4.64 

299.60 

58.80 

9403310413 

-16.08 

6.24 

310.80 

58.80 

9403310517 

-17.52  . 

5.92 

309.40 

58.80 

9403310617 

-18.06 

6.08 

310.80 

60.20 

9403310719 

-19.56 

5.44 

315.00 

60.20 

9403310822 

-19.50 

7.20 

310.80 

60.20 

9403310922 

-20.34 

7.20 

316.40 

60.20 

9403311025 

-21.84 

5.60 

315.00 

60.20 

9403311125 

-22.68 

4.80 

309.40 

60.20 

9403311228 

-23.46 

3.68 

299.60 

60.20 

9403311331 

-23.82 

4.32 

310.80 

60.20 

9403311431 

-24.48 

3.20 

305.20 

60.20 

9403311534 

-25.26 

2.56 

302.40 

60.20 

9403311634 

-24.96 

2.72 

299.60 

60.20 

9403311736 

-24.72 

2.56 

303.80 

58.80 

9403311841 

-23.88 

3.20 

310.80 

58.80 

9403311941 

-23.52 

2.88 

303.80 

58.80 

9403312044 

-22.98 

3.04 

295.40 

58.80 

9403312144 

-22.44 

2.24 

298.20 

57.40 

9403312246 

-21.84 

2.40 

298.20 

57.40 

9403312349 

-21.42 

2.24 

305.20 

57.40 

9404010049 

-20.88 

1.92 

315.00 

57.40 

9404010151 

-20.70 

1.44 

201.40 

57.40 

9404010254 

-21.00 

1.28 

180.40 

57.40 

9404010354 

-21.60 

0.16 

207.00 

58.80 

9404010458 

-22.56 

0.00 

313.60 

58.80 

9404010558 

-22.56 

0.48 

268.80 

58.80 

9404010600 

-22.56 

0.48 

268.80 

58.80 

9404010700 

-23.70 

0.80 

268.80 

60.20 

9404010804 

-24.72 

0.80 

256.20 

60.20 

9404010904 

-24.00 

1.28 

285.60 

60.20 

9404011008 

-24.72 

2.08 

312.20 

60.20 

9404011108 

-24.30 

1.76 

289.80 

60.20 

9404011209 

-25.38 

1.44 

278.60 

58.80 

9404011313 

-24.84 

1.76 

288.40 

58.80 

9404011413 

-24.90 

2.56 

292.60 

60.20 

9404011517 

-25.08 

2.24 

303.80 

60.20 

9404011617 

-24.54 

2.24 

301.00 

60.20 

9404011719 

-24.00 

2.72 

287.00 

58.80 

9404011822 

-23.10 

3.20 

288.40 

58.80 

9404011922 

-22.14 

4.00 

295.40 

58.80 

9404012026 

-20.94 

3.84 

292.60 

58.80 
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YrMoDaHrMi  AirT(C)  WSpd(m/s)  WDir  Comp(Floe  Rotation) 


9404012126 

-19.80 

4.96 

291.20 

58.80 

9404012228 

-19.26 

5.12 

282.80 

58.80 

9404012331 

-19,02 

5.76 

277.20 

58.80 

9404020031 

-18.66 

6.56 

291.20 

58.80 

9404020134 

-18.60 

6.56 

284.20 

58.80 

9404020237 

-18.78 

7.04 

292.60 

58.80 

9404020337 

-19.26 

7.20 

289.80 

58.80 

9404020441 

-20.04 

7.20 

284.20 

58.80 

9404020541 

-20.10 

8.32 

287.00 

58.80 

9404020642 

-20.10 

8.80 

289.80 

58.80 

9404020745 

-19.92 

7.68 

292.60 

58.80 

9404020845 

-19.20 

9.44 

298.20 

57.40 

9404020949 

-19.38 

12.96 

291.20 

57.40 

9404021049 

-19.86 

12.80 

289.80 

58.80 

9404021152 

-19.02 

11.04 

294.00 

58.80 

9404021253 

-18.18 

19.84 

295.40 

57.40 

9404021355 

-18.24 

20.32 

298.20 

57.40 

9404021458 

-17.58 

20.48 

320.60 

57.40 

9404021558 

-17.34 

20.64 

299.60 

57.40 

9404021600 

-17.34 

20.64 

299.60 

57.40 

9404021700 

-16.92 

20.16 

299.60 

56.00 

9404021804 

-16.38 

20.48 

323.40 

56.00 

9404021904 

-15.48 

21.92 

299.60 

57.40 

9404022006 

-14.70 

21.92 

298.20 

57.40 

9404022107 

-13.92 

23.68 

298.20 

56.00 

9404022209 

-13.56 

25.28 

302.40 

57.40 

9404022312 

-13.32 

22.56 

326.20 

57.40 
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